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Abstract 

We study a classical spin model (more precisely a class of models) 
with O(N) symmetry that can be viewed as a simplified D dimen- 
sional lattice model. It is equivalent to a non-translationinvariant one 
dimensional model and contains the dimensionality D as a param- 
eter that need not be an integer. The critical dimension turns out 
to be 2, just as in the usual translation invariant models. We study 
the phase structure, critical phenomena and spontaneous symmetry 
breaking. Furthermore we compute the perturbation expansion to low 
order with various boundary conditions. In our simplified models a 
number of questions can be answered that remain controversial in the 
translation invariant models, such as the asymptoticity of the per- 
turbation expansion and the role of super-instantons. We find that 
perturbation theory produces the right asymptotic expansion in di- 
mension D < 2 only with special boundary conditions. Finally the 
model allows a test of the percolation ideas of Patrascioiu and Seiler. 



1 This work is based to a large extent on the second named author's doctoral 
dissertation (nj 
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1 Introduction 



In their well-known paper on the Mermin- Wagner theorem |l|], Dobrushin and 
Shlosman considered in a side remark a model that is the prototype of the 
model studied in this paper. Their purpose was to illustrate the importance 
of the smoothness of the interaction for the question of symmetry breaking. 
Here we take advantage of the fact that this model on the one hand has 
almost the simplicity of a one-dimensional model but on the other hand has 
a tunable parameter D playing the role of dimension. For integer values of 
D, the model can actually be implemented as a model on the lattice Z D with 
nontranslational Hamiltonian. 

Such a model was also studied in in order to verify that the smoothness 
restrictions of ]l| can be relaxed and are only needed in a neighborhood of 
the minima of the Hamiltonion. Finally, in Georgii's book || a similar non- 
translation invariant chain, but with Ising spins is analyzed. 

Here we use this type of model as a laboratory to test various ideas pro- 
posed by Patrascioiu and the first named author in their quest to prove that 
the conventional distinction between abelian and nonabelian is unjustified. 
The plan of the paper is as follows: After giving the definitions of the model 
and the various boundary conditions used, in Section 3 we use the full, non- 
perturbative solution of the model to study the phase structure as a function 
of the 'dimension' D. It turns out that the critical dimension is still D = 2; 
for D < 2 the model does not show spontaneous magnetization or phase co- 
existence, whereas for D > 2 is does. For D = 2 the model (without external 
magnetic field) and does not show asymptotic freedom. 

In Section 4 we contrast those nonperturbative results with the results of 
(low order) perturbation theory (PT). We find that for zero magnetic field 
PT at the level of one loop already becomes dependent on the boundary 
conditions (b.c.) used for all D < 2 and therefore in general does not yield 
the correct asymptotic expansion of the model. The analogous result for the 
ID model is well known (@, |5|). In this simple class of models it is, however, 
easy to find b.c. in which local observables are independent of the size L of 
the system and hence PT of the finite system does give the correct asymptotic 
expansion. Such b.c. exist in principle also for the full translation invariant 
models. Formally they arise by integrating out all the variables outside a box 
of suitable size, more precisely these b.c. are obtained using the Dobrushin- 
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Lanford-Ruelle (DLR) equations @, 0] . In our casefor 1 < D < 2 those DLR 
boundary conditions turns out to be simply the standard free b.c.. Whereas 
the nonperturbative results are rather easily obtained for our model, PT 
turns out to be harder to compute; for this reason we limit ourselves to one 
loop. This is sufficient, however, for seeing all the effects and subtleties we 
are interested in. 

Finally, in Section 5, we discuss some percolation properties of our model. 
In particular we test the ideas of || Q on the percolation properties of various 
sets defined by spins pointing in certain subsets of the spheres SW-i- We find 
that in 1 < D < 2 we are generally in the situation of 'critical percolation', 
as suggested by [§, [|. 

We should stress that our analysis shows that there is no qualitative 
difference between the abelian case N = 2 and the nonabelian one (N > 2). 
In that sense it lends support to the 'heretical' scenario of Patrascioiu and 
Seiler that predicts the existence of a 'soft' phase in all 2D O(N) models. 
Sceptics might still argue, however, that the model is more ordered than 
the standard translation invariant ones and that this is the reason for the 
existence of the soft phase. 



In this paper we are to a large extent presenting results of [[IIJ, with 
emphasis on the physical interpretation. For more mathematical and com- 
putational details we refer the reader to that work. 



2 The Model 
2.1 Definition 

To motivate the model, we start with a general class of classical O(N) spin 
models defined on lattices Z D , but with link dependent couplings. To each 
lattice site x G TL D we associate a classical O(N) spin, i.e. 

x^s(x), s(x)eR N , ||«(ar)|| = l (2.1) 

The Hamiltonian is of the form 

H = -Y,Mx)-<y)-hY / s N (x) (2.2) 

(xy) % 
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where sjv(i) = s(j) • and e^v is the unit vector in R pointing in the iVth 
direction; the sum is over nearest neighbors and h represents an external 
magnetic field; the dot denotes the standard euclidean scalar product in ~R N . 

We now choose an origin and surround it with a family of concentric 
quadratic (cubic, hypercubic) 'shells' (see Fig.[TJ). We 'freeze' all the links 
(nearest neighbor pairs) sitting inside one of the shells by sending the corre- 
sponding (3 xy — > oo, thereby forcing all spins within such a shell to be equal. 
All other couplings (3 xy are set equal to a constant (3. 




Figure 1: Scheme of the lattice for D = 2; on the thick lines we send (3 xy — > oo 

We can therefore identify all the points within a 'shell'; thus the resulting 
model can be equivalently described as a non translation invariant semi- 
infinite spin chain formally defined by the Hamiltonian 

oo 

H = ~ £(Mi) ■ <J + 1) - h jSN (j)) (2.3) 

3=1 

where 

bj = 2D(2j - 1) D ~\ j = 1,2,3,... (2.4) 
hj = hgj (2.5) 

with 

gj = [(2j-l) D -(2j-3)% j = 2,3,4,... (2.6) 
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and 

9i = I- (2-7) 

Note that for j — > oo both 6j and gj behave asymptotically like 2D(2j) D ~ 1 . 

As usual, we will have to study first finite chains of length L and then 
take the thermodynamic limit L — > oo. This will require imposing certain 
b.c. at the end of the chain. In addition we will also generalize the model 
by introducing similar b.c. at the beginning of the chain. The form of the 
Hamiltonian given above corresponds to free b.c. at the beginning of the 
chain, which is natural from the point of view of the .D-dimensional lattice, 
but we will also be interested in posing Dirichlet b.c. at the origin, which 
in combination with Dirichlet b.c. at j — L will correspond to the superin- 



stanton b.c. introduced in [[TT|| . A class of b.c. that allows to interpolate 
between free b.c. and Dirichlet b.c. is given by the following finite volume 
Gibbs measures: 

1 



dfi L = — e"^ n (2-8) 

3=1 



Z L 



where dv is the O(N) invariant probablity measure on SV-i and 

L-l L 

Hl = ~J2 b Aj) ■ s{j + 1) - E hjS N (j)) - as N {l)) - bs N {L)) (2.9) 

3=1 3=1 

Free b.c. at j = 1 (j = L) are given by putting a = (6 = 0); Dirichlet b.c. 
at j = 1 (j = L) are obtained by taking the limit a — > oo (6 — ► oo). Putting 
b = bi, however, also corresponds to Dirichlet b.c, but at j = L + 1 for a 
chain of length L+ 1. Likewise a = 1 can be interpreted as Dirichlet b.c. on 
an extended chain. 

We will also make use of more general b.c. at j = L which are obtained 
by replacing the factor exp(/3bsisr(L))di>(s(L)) in the Gibbs measure for a 
finite chain with a general positive measure written formally (in distributional 
notation) as ^(s(L))du(s(L)). We may even let ip depend on the length L 
of the chain; in particular this is necessary if we want to describe Dirichlet 
b.c. at j = L + 1 by choosing b = b L , as described above. 

The model we have defined contains the dimension D as a parameter 
which no longer has to be an integer (it could even be chosen complex). In 
the following we will treat D as a real parameter > 1 and we will use almost 
exclusively the representation of the model on a chain. 
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Let us note one crucial fact which turns out to be responsible for the 
dependence of the properties of the model on the 'quasidimension' D: 

Proposition 2.1: For 1 < D < 2 and any k > X/^fc^J 1 — * 00 f° r 
L — ► oo, whereas for D > 2 X/yLfc&J 1 converges to a finite limit. More 
precisely, for D < 2 EjU 6 / 1 = 0(fc 2 - D ), for D = 2 it is O(lnJfe), and for 
D > 2 it approaches its limit like k 2 ~ D . 

The proof is a straightforward consequence of the asymptotics of the bj. 



2.2 Thermodynamic quantities 

To define thermodynamic quantities like free energy density, magnetization 
or susceptibility, we first note that the volume of the chain of length L is 
naturally defined as 

V L = (2L + 1) D . (2.10) 
The free energy density then becomes 

1l = ~^Z l (2.11) 



and the magnetization we define as usual as 

dfi 
dh 



M L = - d 4± (2.12) 



Thus 



1 L 



Ml = t F Y.9 3 (s n {3)). (2.13) 

3=1 

Likewise we define the (longitudinal) susceptibility as 



which means 



L 



Xl = t^J2 9j9l 

Vl 3,1=1 



(snU)sn(I)) ~ (snU)) (sjv(0) 



(2.15) 



All these definitions are analogous to those in translation invariant systems. 
Of course we are mostly interested in the thermodynamic limit L —>■ oo. We 
use the definitions ( |2.13| ) and ( |2.15|) to define the spontaneous magnetization 
and the susceptibility also for the thermodynamic limit at h = 0. 
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2.3 Superinstantons 



For D < 2 the model has superinstantons, just like the translation invariant 



model [pTin - These are configurations of arbitrarily low energy which are 



disordering the system and are responsible for the absence of spontaneous 



symmetry breaking (Mermin- Wagner theorem [|12|1). 

More concretely, these configurations are obtained by imposing b.c. at 
j = i and j = k such that s(i) ■ s(k) = arccos(a) < 1 and minimizing the 
energy under those b.c. The existence of such a minimizing configuration is 
obvious, because the energy is a continuous function on a compact set. It is 
clear that in this minimal configuration the spin will vary on a great circle 
in Sn-i, so we can describe the configuration by the angle 

4>j = arccos (s(j) ■ s(j + 1)) . (2-16) 

The quantitiy to be minimized is then 

E(i,k) = -J2bAcos^ J )-l) (2.17) 

j=i 

under the condition that 

fc-i 



$>i = a (2.18) 

Using a Lagrange multiplicator A one obtains therefore the equations 

A - bj sin(0j) = (2.19) 



which has the solutions 



bj = arcsin(^-) (2.20) 

"3 



and yields for the minimizing configuration 



fe-i 



E sA .(i,k) = J2 b 3 1 



with 

fc-i 



\ 




(2.21) 



arcsin(^-) = a. (2.22) 

j=i h 3 
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We will now derive upper and lower bounds for E S 4_(i, k). Using 

1 - x < VT^x < 1 - I for < x < 1 (2.23) 

we obtain 

i fc-l \ 2 fc-l \2 

£ (2 . 2 4) 

3=1 J 3=i J 

Bounds for the Lagrange multiplier A are obtained using 

x < arcsin(x) < 2x for < x < 1 (2.25) 
from which we obtain, using eq.( p.l9|) 

fc-l k-l 

AE 6 / 1 < a < 2A E 6 j t1 ( 2 - 26 ) 

j=i j=i 

or 

r— < A < —T-, — r . (2.27) 

Combining this with eq.( |2.24|) we finally obtain 

o? „ , , a 2 



sv fc-i,-i ^ fc ) ^ ^FT1-T (2-28) 

Here we can see again the distinction between low (D < 2) and high (D > 2) 

quasidimension: 

For D < 2 

lim E sA .(i,k) = (2.29) 

fc^oo 

whereas for D > 2 E si (i, k) is uniformly in /c bounded away from 0: 

a 2 

E s . L (i,k)> , >0. (2.30) 

The fact that the superinstantons in D < 2 cost arbitrarily little energy is 
responsible for the fact that the system has no long range order (just as the 
translation invariant models are according to the Mermin- Wagner theorem 



T2|). This will be proven in the next section. 

Due to the freezing of all the spins within one 'layer', the model does not 
have instanton-like configurations. 
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3 Nonperturbative Analysis 



3.1 Transfer Matrices 

Because our model only couples nearest neighbors, it can be described easily 
in terms of transfer operators, which are of course site dependent. These 
transfer 'matrices' are trace class operators on the Hilbert space 

H = L 2 (S N - U dv), (3.1) 

which we normalize such that the largest eigenvalue equals 1. In slight abuse 
of notation, we will use the same symbol for the operators and their integral 
kernels. 



The normalized transfer matrix from site j to site j + 1 is 

(3(b jS (j) ■ s{j + 1) + h jSN (J)) 



2}(*0'),*0' + l)) = -exp 



Z 3 



(3.2) 



with Zj chosen such that 



l?ll = l (3-3) 



These transfer operators are not self-adjoint, but it is easy to transform them 
by a similarity transformation with a bounded operator into the self-adjoint 
operators with integral kernel 



1 



h 



Ti(8(j),8(j + 1)) = -exp (3(b jS (j)-s(j + l) + ^(s N (j)+s N (j + l)) 



(3.4) 



These operators are in fact positive, as can be seen easily be expanding 
exp[s(j) • s(j + 1)]. 

Since the transfer operators have also positive integral kernels, by a trivial 
change of normalization it is possible to interprete them as transition prob- 
abilities and the whole system as a Markov chain. This point of view will, 
however, not play a great role in this paper. 

We introduce a shorthand for products of transfer matrices: 

k— 1 

T jk = J] % (3.5) 
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and define 



(3.6) 



Using the transfer matrices, it is possible to rewrite for instance the ex- 
pectation value of a spin component at site k 



(s a (k)) L = / djJt, L s a {i) (3.7) 



in the form 



(sMl = ^ TuSa%L2p) (3.8) 



where we denote by (.,.) the scalar product in TC. ip and ip are suitable 
vectors in TC = L 2 (Stv-i) describing the b.c. (actually we may even replace 
them by general positive measures on SV-i)- The spin variable s a appearing 
on the right hand side of eq.( |3.8|) is to be interpreted as a multiplication 
operator on L 2 (Sn-i)- 

The spin-spin two point function 



(s a ({)s b (k)) L = I dfi L s a (i)s b (k) (3.9) 
for i < k can similarly be expressed as 

(s a (t)s b {k)) L = (3.10) 

For a nonvanishing magnetic field we do not expect any interesting phe- 
nomena; the transfer matrices will force the spins 'at infinity' to be aligned 
with the direction of the magnetic field e^, therefore one expects a unique 
thermodynamic limit, independent of the b.c. imposed at L. 

For vanishing magnetic field the situation is more interesting; therefore 
from now on we will assume h = 0. To analyze the situation, we need some 
preparation (see for instance |13], [TJJ): 

The Hilbert space TC = L 2 (Sn_i, dv) can be decomposed into the eigenspaces 
TCi of the Laplace-Beltrami operator Alb on the sphere S^-i: 

oo 

H = ®Tb (3.11) 
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The projections V\ onto the eigenspaces are integral operators which for 
N > 2 can be expressed in terms of the Gegenbauer polynomials Ci (see 



15,10 



n , i H r (^_i + /_ m 

C? \x) = —t- t rT(-l) m V ,„ „ \ i y (2x)'~ 2m , iV>2, 

1 1 ^ r(f -l)if ^ m!(/-2m)! 1 ; 

^ VU — m) 

C?(x) = V(-ir^7 , ' r(2x) z - 2m , 1^0. (3.12) 

iV ; r(m + l)r(/-2m + l) V ; T K J 



as follows: 



21 -\- N — 2 v i 
Vi{a,J) = M+ N _ 2 Cr\s-s') (3.13) 



For N = 2 the Gegenbauer polynomials have to be replaced by the Chebyshev 
polynomials of the first kind: 

T (x) = 1 (3.14) 



1 



N 



1 , 



TAx) = lim — C, 2 (x) fori > 1. (3.15) 

^ ' N ^ 2 N — 2 y J 

In the following we will write the equations in the form valid for N > 2, 
involving the Gegenbauer polynomials, with the understanding that for N = 
2 the analogous formulae involving the Chebyshev polynomials hold. 

Due to the O(N) invariance of the transfer matrices for h = 0, all Tj 
commute with those projections, and the spaces Tii are also simultaneous 
eigenspaces of all Tj. Because the Gegenbauer (Chebyshev) polynomials form 
a complete orthogonal set on the interval [—1, 1], the integral kernel of the 
transfer matrix can be expanded in the sense of L 2 ([— 1, 1], (1 — x 2 )^~dx) as 

oo 

T 3 {s ) s') = Y,c l C? /2 -\s.s') (3.16) 

1=0 

From this we see that each subspace Hi is an eigenspace with eigenvalue A/ 
of Tj. The eigenvalues A; are therefore given by 

Xi(Pi) = — h^ = — dxexp(/3,x)(l -x 2 ) 2 /, , 
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1 N-l /2\t* 

where J r (.) is the modified Bessel function ( ||15| ). From the fact that the Tj 
are positivity improving operators (see |17| ) it follows that 



< Xi(Pj) < X (Pj) = 1. (3.18) 

The normalization X (/3j) = 1 allows to determine Zj explicitly and we get 
for the eigenvalues finally 

m) = j 2 (n \ ( 3 - 19 ) 



3.2 Thermodynamic limit 

In this subsection we will discuss explicitly the thermodynamic limit for the 
one and two point functions of the spins. The generalization to higher n-point 
functions is straightforward in principle. 

To take the thermodynamic limit we need some information about the 
limit of the product of a large number of transfer matrices. Since for h = 
all transfer matrices commute and have the same eigenprojections Vi, we 
obtain 

oo k— 1 

T ik = E II (3.20) 

1=0 j=i 

So we have to study the behavior of the products YljZi Xi(f3j) for large k. 



For this purpose we can use some results of []14| 1131 , where we analyzed the 
asymptotic behavior of such eigenvalues (actually for more general forms of 
the Hamiltonian) . The main result is 

In m) = 1 - Kl + "~ 2) + 0{?f) (3.21) 

We define for i < k 

k—l 

At(i,k)= UWi) (3-22) 

j=i 
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Using ( |3.21| ) and Prop. 2.1 we find the asymptotic behavior of In A;(l, k): 
v ' ; 2(3 2 D D(2-D) \ (3 2 



^ + ^-2)1 , „ /a _ 2 , 



(3.23) 



lnA z (l,fc) = — J --\nk + 0{(3~ 2 ) for D = 2 (3.24) 



and 



in A,<», oo) = - ' ( ' + * - 2) ^ _ 2) + O (*^) fo rfl >, 

(3.25) 

Using this, we obtain easily 

Proposition 3.1: For D < 2 = s-limf.- >oa Tii t = "Po- For D > 2 7^ 
exists but is strictly positive; s-lim^oo 7^ = t. 

Proof: To prove strong convergence, it suffices by a 3e argument to prove 
convergence of all the eigenvalues. We now treat three cases separately: 

(1) D = 1: All Tj are equal and the result follows trivally from the fact that 
for I > 1 A* < 1 

(2) 1 < D < 2: The assertion follows from the fact that for all / > 1 

lim A,(l, k) = 0. (3.26) 

k^oo 

To see this, use 

Inll Wi) < E(W) - !) = -E[ ^ + 2 ^~ 2) +Q(^ 2 )]- ( 3 - 27 ) 

j=i j=i j=i Pj 

According to Prop. 2.1 the first sum diverges to — oo; more precisely, it 
behaves like — (k) 2 ~ D for 1 < D < 2 and like In k for D = 2. The sum of the 
correction terms diverges for 1 < D < 3/2 at most like k 3 ~ 2D , for D = 3/2 
like In k and converges for D > 3/2. In any case the asymptotic behavior for 
large k is determined by the leading term; so the product of the eigenvalues 
diverges to zero and the assertion follows. 

(3) D > 2: The product of the eigenvalues Ai(k, oo) converges absolutely to 
a nonzero value because 

k 

EKWi)- 1)1 (3-28) 

j=i 
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converges for k — > oo. ■ 

One and two point functions simplify in the thermodynamic limit, pro- 
vided either 1 < D < 2 or we have free b.c. at 1 (i.e. a = 1. Free b.c. at 
1 seem most natural anyway from the point of view of the D dimensional 
lattice; if we consider O(N) invariant observables, however, it is just as le- 
gitimate to choose a ^ 0, since in the thermodynamic limit this only means 
fixing the global O(N) invariance. 

We first discuss the case 1 < D < 2 and a = 0. In this case the Hilbert 
space vector i/i in the expression ( [3781) becomes proportional to the 'ground 
state' ipo, which is the function identically equal to 1 on Sjv-i- i>o spans the 
range of Vo', it is invariant under all Tj and therefore we obtain 

Mi)h = ^TTT 1 < 3 ' 29 > 

Wo, W) 

and 

(s a {t)s b {k)) L = — T (3.30) 

Next we use the fact that s a maps the subspace Ho into the subspace Hi, to 
conclude (assuming i < k) 

(s a (i)) L = ^0-Mi,L) (3.31) 

and 

{s a (i)s b (k)) L = — ~- Ai{t, k). (3.32) 

Now we take the thermodynamic limit, using Proposition 3.1. We obtain for 
1 < D < 2: 

(aa(«)>oo = ( 3 - 33 ) 

and 

(s a (0 s 6(^))oo = (V'o, s a s 6 V>o)Ai(i, k) (3.34) 
By O(N) invariance it is easy to see that 

(ipo,s a s b ip ) = ^S ab (3.35) 
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More generally, by the same reasoning, if F is a bounded measurable function 
of finitely many spins s(ii), ...,s(ik), and F its O(N) average, 

^>oo = (^> 00 - (3-36) 

Finally let us discuss what happens for 1 < D < 2 and a ^ 0, provided 
we consider ON) invariant observables. So let F now be an O(N) invari- 
ant bounded measurable function of finitely many spins. In this case the 
Hilbert space vector ip in the expression (|3.8|) becomes a function ip{ a ) with 
{ip{a), ipo) 7^ 0. But in the thermodynamic limit, again using Propositionn 
3.1, we obtain 

-v (il>(a),Fip ) - 

F ) i = < \ = wo> F w, 3.37 



OC.i 



which is independent of a and 6. 

We summarize what we have found for 'low' dimension in 

Theorem 3.2: For 1 < D < 2 and a = (free b.c. at 1) in the thermo- 
dynamic limit, irrespective of the b.c. at L (i.e. b) the following holds: 

• (*.(*)>«, = o 

• (*«(*)*&(*;))«, = ^<UIl*=i AiOfy) = jjSabA^k). 

In addition for any bounded measurable function F of finitely many spins, 
and F its O(iV) average, we have 

(^)oc = (f)„ ■ (3.38) 



Furthermore for arbitrary a (F) is independent of a and b. Finally for 

\ / oo 

k — > oo we have the asymptotic behavior 

<«.(0*(AO>oe ~ ^ ex P (- 2 D2^(2-D) fc2 " C ) (1 " D < 2) (3 - 39) 
and 

<aa(0^(*)>oo~^«pf-^llnA;J (£> = 2). (3.40) 
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For D > 2 the situation is different and a little more involved. We have 

(3.41) 



<*a(*))oo = -T, 7-MhOO) 



and 



(s a (i)s b (k)) c 



(jj , s a s b T koo jj) 



Ai(i,fc). 



(3.42) 



If we choose Dirichlet b.c. by sending ip to a delta function concentrated 
at s = ejv, the one point function simplifies to 



<«.(<)> 



oo,Dir 



5 aJV Ai(i, oo) 



(3.43) 



We can analyze the 2 point function further by using the fact that only 
the subspaces Ti Q and 7i 2 contribute here; generally we have 



Sab 



N 



(3.44) 



for any G 7i. Inserting <p = T^oo^ we obtain 



(saitjsbik)}^ = o ab —— — ~—Ai{%, k)+ — — Ai(z, k)A 2 {k, oo). 



N (M) 
Using the fact that 

(V'O,'!/') 

and taking the limit fc^oowe obtain 



lim (s a (i)s b (k)) 



fc^oo 



1 1 



(3.45) 

(3.46) 

Ai(i,oo). (3.47) 



This latter expression is in general nonzero, but for Dirichlet b.c. (tp = 5e N 
we can evaluate it further, because in this case 



(ipo, (s a s b - -^s 2 )^)) = 5 aN 5 bl - 5 ab 



and we obtain 



N 



lim (s a (i)s b (k)) Dir 

fc— >oo 



1 

TV 



5 a N$bNAi(i, OO) 



(3.48) 



(3.49) 
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Comparing with the result ( |3.43j ) for the one point function we see that the 
latter limit equals 

lim (saii))^ Mk))^ (3.50) 

fe— >oo ' ' 

In other words, for Dirichlet b.c. at L — > oo, the truncated two point function 

(S a (l); Sbik))^ ^ = (Sa^Sbik))^ ^ - (s a («))oo,Dir ( S b( k )) oo,Dir ( 3 - 51 ) 

goes to for k — > oo. In this sense we may interprete the state with Dirichlet 
b.c. at oo as a pure phase of the system. Note also that by fl3.25| ) the limit 
is approached like k 2 ~ D . 

We can still work out the truncated two point function for Dirichlet b.c. 
at oo in a little more detail: from eq. ( |3.45|) we have 

(s a (i); s 6 (A;)) 00)Dir = 5 ab -^Ai(i, k) + (s aN 8 bN - 8ab-j-) Ai(i, k)A 2 (k, oo) 
- 5 aN 5 bN A 1 (i,oc)A 1 (k,oo) (3.52) 

which can be rewritten as 

SaN$b N ~ da,bj-f) Ai(i, k)(A 2 (k, oo)-l)+5 aN 5 bN (A x (i, k) - A x (i, oo)Ai(/c, oo)) 

(3.53) 

For a = b ^ N this is manifestly positive. Furthermore, it is also straightfor- 
ward to see that the O(N) invariant truncated two point function is positive: 

(s{i);s(k)) > 0. (3.54) 

where the semicolon is meant to include the dot symbolizing the scalar prod- 
uct in R N . 

Let us now summarize what has been found for D > 2 in 

Theorem 3.3: For D > 2 and and a = (free b.c. at 1) in the thermo- 
dynamic limit 

• the one-point function is 

(SaW)oo = Ai(«,oo)— — ~— (3.55) 

which is nonzero provided the scalar product in the numerator does not 
vanish. This is in particular the case if s a = stv and the b.c. parameter 
b > 0. 
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Hindoo (s a (i)s 6 (A;)) oo is generically nonzero. 

For Dirichlet b.c. (b — > oo) the one point function is 

( S ^))oc,Dir = $aNM(i, Oo); (3.56) 

furthermore the two point function has the cluster property 

lim («a(i);* 6 (*)> oo ^r = 0. (3.57) 

k — >oo 



The truncated two point function behaves like 

^ oo,Dir 



( Sa (i); Sa (k)) EE M»>6(*)>ooJ» - MOJocDi. ( S »(^))oc,Dir = O(^) 



(3.58) 

From the two theorems we can easily derive results about the magnetiza- 
tion: 

Cor. 3.4: For 1 < D < 2 the spontaneous magnetization = 
lim^oo Ml vanishes, irrespective of the b.c. For D > 2 the model has 
spontaneous magnetization with suitable b.c; for Dirichlet b.c. = 1. 

Proof: The spontaneous magnetization is, according to |2.13| a certain 



weighted average of the 1-point function (sjv(«))oo- The one P om t function 
converges for % — > oo to 

= %^ (3.59) 

Any average of a convergent sequence is equal to its limit, so we have 

Moo = (sjvM)oo (3- 60 ) 

For Dirichlet b.c, because ip becomes the distribution 5e N concentrated at 
ejy, the last quantity is equal to 1. ■ 

The physical interpretation of the fact that the truncated two point func- 
tion (s a (i)s a (k)) decays faster than any power of k for 1 < D < 2 is that 
the model is in its high temperature phase for any value of f3. Likewise 
the power-like decay in D = 2 means that the model is critical for any (3. 
The power-like decay of (s a (i); s a (k)) for D > 2 expresses the presence of a 
Goldstone-like mode in the system. 
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The susceptibility is a little pathological in our model. As defined above, 
it diverges in general (except for D = 1) due to the freezing of the spins 
within one layer. Namely, if the truncated two point function is nonnegative 
(as it is, according to the discussion above), we have 

XL = jriZ 9i9k (s N (i); s N (k)) > £ g? (s N (i) 2 ) = -r^r X> 4 2 (3.61) 
which diverges for L —>■ oo. 



3.3 Presence and absence of asymptotic freedom 

In this subsection we discuss the issue of asymptotic freedom using a defini- 
tion of the Callan-Symanzik /^-function Pes due to Patrascioiu ||. We use 
the truncated two point function 

G p [i,k) = (s(i);s(k)) (3.62) 

to define the following Renormalization Group (RG) invariant quantity: 

F ^ k) = ^kJ (3 - 63) 

We rescale the lattice distances in Fp and ask how this can be compensated 
by a change of the coupling constant g = 1/ y 7 /?. 

This compensation cannot be made exact by only changing (3, but it 
works aymptotically in the limit of many iterations, or equivalently for the 
large distance asymptotics. 

Since we have the asymptotic behavior (see ( |3.25| )) 



(«« • - ex P ( 2 dd(D-2)0 " k2 ~°)) (3 - 64) 

(for k > i), it is seen straightforwardly that a doubling of i and k can be 
compensated by replacing (5 with 2 2 ~ D ' (3 or g with 2 D I 2 ~ 1 . Of course there are 
corrections of order 1/ f3 k to the exponent and we will compute the first one 
of these corrections in the next section, but they do not affect the qualitative 
conclusions. 
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If we interpolate Fp to obtain a smooth function on R 2 we can do this 
rescaling infinitesimally and find that Fp obeys, at least asympotically, the 
following RG equation: 

(d t + Pcs(g(t))d g ) F p {eH, e%) = (3.65) 

This equation should really be interpreted at the definition of the Callan- 
Symanzik /3-function, i.e. we have to set 

PcsW)) = - d g F p {e%e% (3 - 66) 
which yields, putting t — and g(0) = g, 

Pcs(g) = ^Y^g + 0(g 2 ). (3.67) 

This result shows that for D < 2 the model is asymptotically free, whereas 
for D = 2 it is critical for any f3 > 0, i.e. we have a half-line of fixed points. 
This is true for any N > 1, i.e. there is no qualitative difference between the 
abelian and the nonabelian versions of the model. 

For D > 2 the analysis is a little different, because we have to take the 
truncation into account. The asymptotic behavior of the truncated two point 
function is 



W);*W)~- V p~- 2) / - D < 3 - 69 > 

Forming the renormalization group invariant Fp we see that the dependence 
on (3 as well as the scale parameter drops out. The RG equation is satisfied 
with 

Pes = (3.70) 
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Of course, strictly speaking, j3 C s is left undetermined by eq. ( |3.65| ) and eq. 
( |3.66| ). But Pes — also is the right answer if we consider the more general 
RG equation satisfied by the truncated two point function itself: 

(d t + Pcs(g(t))d g — (D — 2)) G p (eH, e%) = (3.71) 

We conclude that the model is thus also critical for D > 2; the reason is of 
course the presence of the Goldstone-like mode. 

3.4 Gibbs states and phase structure 

In this subsection we will discuss the set of Gibbs states of our model for 
the semi-infinite chain, obtained by taking the thermodynamic limit of finite 
chains with various b.c. As before, we will use free b.c. at 1 (a = 1), but 
the following should be noted: 

Remark: If a Gibbs state (.) is O(N) invariant, we can modify it by 
introducing an arbitrary measure for a particular spin, without affecting 
the expectation values of O(N) invariant observables. In particular these 
expectation values become independent of the measure ip chosen for the b.c. 
at 1. 

Gibbs (= equilibrium) states can be characterized by the DLR equations 
H |7j. Consider a local observable 0(s(i), ...s(k)), i.e. a continuous function 
of a finite number of spins (1 < i < k); then the DLR equations imply that 
for any integer r > 

(0) = (6 r ) (3.72) 
where O r is a function of s{i) and s{k + r) 

O r (s, s') = ^J 0(s(i), ...8(h)) *n 1 T 3 (s(j), s(j + 1) "llM*®) ( 3 - 73 ) 

j=i j=i 

with 

„ fc+r-l fc+r-1 

Z = j II T s (a{j),a{j + \) II du{a{j)) (3.74) 

j=i j=i 

So among other things, going from O to O r attaches a product of r transfer 
matrices to the observable. We can now ask what happens if we send r —>■ oo. 
In low dimensions (1 < D < 2) that product converges in the strong topology 
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to Vo, as we saw. This is already sufficient (by letting the adjoint of those 
operators act 'to the left') to conclude that for any Gibbs state and 1 < D < 2 

(O) = (OV ) (3.75) 

where we use same symbol P a l so f° r the integral kernel of the operator, i.e. 

V (s,s')=M^)Ms') = l (3.76) 

The insertion of P has the same effect as free b.c. at a point m with m > k, 
so that we conclude 

(0) = (0) iree . (3.77) 

So remembering that we made the general assumption a = 1 in this 
section, we have obtained the result: 

Theorem 3.5: For 1 < D < 2 with free b.c. at 1, there is a unique Gibbs 
state obtained as a thermodynamic limit. In particular this Gibbs state is 
O(N) invariant. 

Remark: As noted above, this implies that for O(N) invariant observables, 
all b.c. at i — 1 are equivalent. 

Now we turn to the cae D > 2. This turns out to be a little trickier 
and we do not obtain a complete rigorous description of all the Gibbs states. 
There is, however, a very natural conjecture. 

Let us first state what we can prove: 

Theorem 3.6: For D > 2 there is an uncountably infinite set of Gibbs 
states parametrized by the set of probability measures on the sphere Sjv-i- 
The extremal states in this set are given by the probability measures 5 e 
concentrated in one point of e e Sn-i- 

Proof: It follows from Proposition 3.1, as in the proof of Theorem 3.3, that 
the thermodynamic limit with a boundary measure ip exists for any local 
observable. Theorem 3.3 says first of all that the one point function in 
general will be different for different choices of ip. We can easily generalize 
the Dirichlet b.c. discussed there by choosing 

i> = 5 e (3.78) 

where Se is the delta function concentrated on a general point e e Sjv-i- 
The one point function will then be 

(s«>oc, e = Ai(*,oc) e (3.79) 
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so all these states are different for different choices of e. 

It is also easy to see that these generalized Dirichlet states are all extremal 
in the set of states given by boundary measures ip, because they satisfy 

(«(«0 ■ e >co,e = 1 (3-80) 
whereas in any other state one has 

(s(oo)-e)<l (3.81) 

■ 

What we do not know rigorously is whether we have exhausted all the 
Gibbs states by our b.c. if), but it is very plausible that this is so. So we 

state 

Conjecture 3.7: For D > 2 the set of Gibbs states is given by the set 
of measures ip on the sphere Sat-i- 



4 Perturbation Theory 
4.1 Preliminaries 

Perturbation Theory (PT) is nothing but the application of Laplace's method 
to the Gibbs measure. For a finite system, the Gibbs factor is very sharply 
peaked around the ground state configuration(s). To make PT work, we need 
a unique such ground state, and we achieve that by a suitable choice of b.c. 
The Gibbs measure is then to lowest order ('tree level') approximated by a 
Gaussian centered at the ground state, and a sequence of corrections to the 
Gaussian arises naturally by Laplace's method. 

For a finite system one can show easily that the resulting expansion in 
inverse powers of j3 is asymptotic to the true expectation values. The usual 
formal PT procedure takes the thermodynamic limit term by term and hopes, 
if that limit exists, to obtain an expansion that is asymptotic to the infinite 
volume Gibbs state. It is well known that this hope fails in some cases (see 
0, H|, [Tl| , and here we will find that it fails in general for dimension D < 2. 



From now on we assume that we have b.c. characterized by b > and 
a > 0. Then the ground state configuration is unique and is described 
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by s(j) = e^v for all j. For large (3 with high probability the spins will 
fluctuate not very far from e^. This motivates the introduction of Cartesian 
coordinates on the sphere Sat-i (as in the classic paper of Brezin and Zinn- 
Justin HIS!) to describe these fluctuations: 



s(j) 



(3- 



2 71% 



with Oj := d=(l — ^ _1 7rJ) 3- 



(4.1) 



These coordinates are singular at the equator, but since we are interested 
in an asymptotic expansion in powers of 1//3, we can ignore this fact. We 
can actually limit ourselves to integrating over the upper hemisphere, i.e. 
always choose the + sign in the definition of Gj. Finally we can extend 
the integration over each irj to all of IR^ -1 . All these changes have only an 
exponentially small effect (in (3) on the integrals, as long as we work with a 
finite system L < oo, and therefore they do not affect PT. 

The partition function can therefore be replaced by 
Z L (h,a,b,D) = 

j ex P [e b 3 (tt, ■ tt, +1 r x 2 ) 1 - p~ l *uf) 

L i i 

+ Ph E 9j (l - 5 + /3a (l - /T 1 ^) 1 +(3b(l- /T 1 ^ 



x exp 



-;I>(i-/r'«J 



7T,-. 



(4.2) 



and accordingly for the expectation values. In other words, for the purpose 
of PT we are reduced to studying the Gibbs measures 



d[i L (ir 



i> •••) 



77 L 



exp(A L ) Yl d-n-j 



with 



(4.3) 



a l = [E h Ui ■ + /? (i - /r 1 *?) 1 (i - /rVJ+i) § ) 
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- ^ 5> (l - ( 4 - 4 ) 

j'=i 

It is now clear how to proceed; we expand A in powers of /3; in this paper we 
will be content to do this to order 

A = 2 D - 1 /3\/ L (l + h) + aPh + b(3 + AW + A (1) + 0(/T 2 ) (4.5) 

where the piece O((3 ) is quadratic in the 7Tj variables: 

^ (0) = -IjLbfa+i - *j) 2 - \J.W) - \*l - \*l (4-6) 



and the term 0{(3 ) is 



A™ = -- L 

2/5 



+ iE^ 2 ) 2 -E^ + i«(^ 2 + ^!) 2 l- (4.7) 

As long as not a = b = h = 0, — A( ) is a positive definite quadratic form in 
the 7T.,- which we write as 

A^ = ~(n,Qn) (4.8) 

where 7r stands for the vector (vrf )"~ 1 1 ''"'^ v_1 in R^ -1 ^ and Q has matrix ele- 
ments Qll = 5 ab qik- We combine A^ with the Legesgue measure to produce 
the Gaussian probability measure 

1 1 

dUcin) = exp(--(7r, Qir) (4.9) 

with covariance C = Q" 1 . C is proportional to the identity operator in the 
internal space R^ -1 , i.e. its matrix elements are of the form Cft = 8 c ik . 

even though it contains also terms quadratic in 7r, is treated as an 
interaction because it is of order j3~ x . 
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4.2 Spin two point function 



We first derive the explicit form of PT up to order f3~ 2 for the invariant 
spin-spin correlation in terms of the covariance of the Gaussian measure up 
to order f3~ 2 . First note that 

(s(i) ■ s{k)) = 1 - JL ((^ - n k f) - _L_ ((^ - ttQ 2 ) + 0(/T 2 ) (4.10) 

provided we have a,b,h > and at least one of them > 0, because then 
the ground state will have all spins aligned parallel to e^. As usual, PT 
is generated by expanding also the interaction in inverse powers of j3 in the 
Gibbs measure; as is well known, the correct normalization leads to the 
phenomenon that 'vacuum graphs cancel', i.e. only terms survive in which 
the interaction is connected by 'lines' (covariances) to the observable. Thus 



(«w ■ < k )) = 1 - 4 ((*< - ^) 2 > c - -L (0* 2 - 




-7r k ) 2 ;A^j + 0(/T 3 ) (4.11) 

where the semicolon indicates that only contributions are to be taken that 
connect the expression to the left with that to the right of it. (.) denotes the 
Gaussian expectation value with covariance c and is defined in eq.([0). 

The term order f3~ x is therefore 

(s(i) ■ s(k)) {1) = —Z±{c« + c kk - 2c ik ) (4.12) 

To order f3~ 2 we find 

(s{{) ■ s(k)) {2) = (/) + (//) + (///) + (IV) + (V) (4.13) 



with 



N-l 



(cl+c 2 kk )(N+l)-Ac 2 k -2(N-l)c ti c kk 

(4.14) 



(II) 



4/3 2 



3=1 



N 



2(3 2 



J^(Qj C k j) 



(4.15) 
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L-l 



N — 1 L_1 f 

EJ(JV + i) 



2 

i+i 



yCjj{Cij 



Ckj) 2 + Cj+l,i+l(Ci,j+l — Cfc,j+l) 2 



" (N-l) 

~~ ^ C j,j+l( C ij ~ C kj){Ci,j+\ — Cfej+l) f , 



(4.16) 



h L 



n2X (jv-iKjv + i) * _ 2 

4^2 £^9j c jj{ c ij c kj) 

(4.17) 

and finally the boundary contribution 



^ ^(tt, - 7r fc ) 2 ; a (tt 2 ) 2 + 6 (tt|) 

(AT — 1)(-/V + 1 



8/5 2 



acn(cn - Cfci) + bc LL (c iL - c kL ) (4.18) 



It is advantageous to rewrite contribution (///) as follows: 

(///) = (///a) + {I I lb) + (///c) 

with 



(4.19) 



(TV — l)jV L_1 r 

(4.20) 

N — 1 L ~ x 

( IIIb ) = —^T h A c n + c i+i,i+i) (cij - °kj - Cij+i + Cfcj+i) 2 (4.21) 



4/3 2 



iV - 1 L_1 

( ///c ) = Yl b j( c jj + c j+iJ+i- 2c jJ+i)( c ij- c kj)((c itj+1 -c kJ+1 ) (4.22) 

i=i 



2(3 2 



4.3 The covariance without magnetic field 

To proceed, we have to compute the covariance C more explicitly. We only do 
this for the simplest case of free b.c. at 1 and no magnetic field, i.e. b = h = 0, 
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but with a > 0. We can ignore the internal space in this computation, i.e. 
put N = 2. Then, with x e R L , we have 

L-l 

(x, Qx) = ]T fyfo+i - x 3 f + ax? (4.23) 
j'=i 

For convenience we define in the following b = a and b L = b. Q is already 
given in the form 

Q = L T L (4.24) 
where L is a lower triangular (actually bidiagonal) matrix with elements 



hk = \Jb k -\ and 

lk+i,k = -yjh for fc = l,2,...,L (4.25) 

all other elements being zero. We now split 

L = D + N (4.26) 
where D = diag(v / &o? •••v^i-i) anc ^ ^ * s nilpotent. Then 

C = Q- 1 = (1 + D- 1 ^)- 1 /}-^! + A^ZT 1 )- 1 (4.27) 

D~ 1 N has only nonzero matrix elements djfc for % — k + 1, and they are all 
equal to —1. Therefore F = (1 + .D^iV) -1 has the matrix elements 

y ik = 1 for i > k 

Uik = otherwise (4.28) 
So, using the shorthand = mm(i, k) we obtain for the covariance 

L m ik -l 

c*fc = E & 7-W% = E bj 1 ( 4 - 29 ) 

j=l r=0 

for 6 = 6^ = 0. If we introduce the further shorthand 

fe— i 

B ik = '£bJ 1 (4.30) 

we can write the covariance we found as 

c ik (a, 0) = B 0m . k = B lm . k + cT 1 (4.31) 
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where the second argument is reserved for the parameter b. 

The covariance for general values of a and b (still without a magnetic 
field) is obtained from this by realizing that changing the b.c. at L is a rank 
one perturbation: 

c ik (a, b) = c ik (a, 0) - (V 1 + c LL (a, 0)) c iL (a, 0)c kL (a, 0). (4.32) 

This can be written in more compact form if we use the definitions M ik = 
max(z, k) and 

B = (a" 1 + r 1 + B.l)- 1 = — !— : (4.33) 

as 

c ik (a,b) = BB 0mih B MihjL+1 (4.34) 
We note some special cases: 

c ik (0, b) = B MikL + ft" 1 = B Mjt , i+1 (4.35) 

c ik (oo,b) = BB lmih B MihjL+1 (4.36) 

c ifc (a, oo) = BB MikL B 0m . k (4.37) 

We also note the following combination of covariances: 

Cjj(a, 6) + c fcfc (a, b) - 2c ik (a, b) = B ik (l - BB ik ) (4.38) 

The behavior for large L and the asymptotic behavior of the infinite 
volume covariances for k — > oo can be obtained from Proposition 2.1, which 
implied that for 1 < D < 2 B ik = 0(k 2 ~ D ), for D = 2 £ ife = 0(\n(k)) and 
for D > 2 Bi k = 0(1). The thermodynamic limit is particularly easy to take 
for b = 0, because in that case by eq. (|4.31|) the covariance does not show 
any dependence on L. 

In general for 1 < D < 2 

lim c ik (a, b) = B Qm (4.39) 

independently of b, provided a > 0. For 1 < D < 2 and a = the thermo- 
dynamic limit of the covariance does not exist. 

For D > 2 the thermodynamic limit always exists, but it depends on the 
b.c. parameter b. 

For the general case with a magnetic field there are no such simple ex- 
pressions for the covariance. A very detailed discussion with many explicit 
and lengthy formulae can be found in |1C]. 
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4.4 Explicit evaluation to one loop 



We will first discuss the simplest case of the 'energy', that is the two point 
function of two neighboring spins. 

To order ('tree level'), we have 
(s(i) ■ s(i + 1)) = 1 - -^{cu + Cj+i.i+1 - 2c iji+1 ) + 0(/T 2 ) (4.40) 



Using eq. (|4.38|) we thus find 



s(i) ■ s(i + 1)) = 1 - ^— - (l - Bbr 1 ) + 0(/T 2 ) (4.41) 



2fl 

and taking the thermodynamic limit of the first order term we obtain 

(s{i) ■ s{i + =br 1 (1<D<2) (4.42) 

and 

(s(i) • s(i + 1)) (1) = bf (l - b Jl Boo ^ j (D > 2) (4.43) 

This means that for D > 2 already at tree level PT shows a dependence on 
b.c. even for such a simple O(N) invariant observable. By expanding the 
nonperturbative solution of the model, analyzed in the previous section, one 
can see that this is a real effect, related to the occurrence of spontaneous 
symmetry breaking (SSB). The dependence on the b.c. parameter b drops 
out, however, if we use free b.c. at 1 (a = 0). This happens because for 
an invariant observable we can always introduce an arbitrary probability 
measure for one and only one spin, without any effect, as explained earlier. 

To the next order ('one loop level') the computation becomes a little more 
involved. Since for D > 2 already the tree level term depends on the b.c, 
we assume from now on 1 < D < 2. 

We first consider the case a > 0, b = 0. It is easy to see that in this case 
term (V) vanishes in the thermodynamic limit. Term (IV) is absent because 
we put the magnetic field h equal to zero. For the other terms, after some 
trivial algebra one obtains 



N - 1 

( / ) = -^^[ Ar + 1 + 4 W, (4.44) 

oft 
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(II) 



JV-1. . 
S-[4(L-<)] 



< ) = " 1) + 4(L - i) + 46^0*] , 



which adds up to 



(s(z)- S (* + l)> (2) (a,0) 



(AT — 1)(N — 3) 
8^ 



(4.45) 
(4.46) 

(4.47) 



The case h = with general b.c. can also be worked out explicitly. First 
note that the term (V) still does not contribute in the thermodynamic limit: 
we have 



00 



(iV-l)(iV+l) 
8p 2 



och(qi - c kl ) 2 + bc LL (c iL - c kL ) 



(4.48) 



For a > and any b the first term goes to as L — > oo, since by eq. (|4.39 ) 
limi-j.oofcji — Cfci) = 0. For the second term we notice 



bc LL (c iL — c kL ) = bB B 0m B m:L+ i(B 0i — B ok ) b = bB B 0m B mtL+1 B ik 







(4.49) 

for L — > oo. 

Inserting our formula eq. (|4.34j) into the expressions at the end of Subsec- 
tion 4.2, the result for (II) and (III) splits into sums from 1 to i — 1 and 
from i to L — 1 or L. For the finite sums the limit L — > oo can be taken 
termwise and is actually zero, because each propagator carries a factor B 
which goes to as L —>■ oo. After some algebra we are left with 



(i) = 

(II) 



N - 1 

-2- [N+1 + 46i5oi] 



8$ 
iV — 1 



b 2 E 



J=i+1 



(Ilia) 



(N — l)N 
4^2 

(I I lb) ~ 



L-l 



67 2 + B 2 Bl £ [^B hL+1 - ABBl 

N - 1 1 

4/3 2 



L+l 



K 2 - ^ E ^fBl L+1 ) 
j=i+l 



(4.50) 
(4.51) 
(4.52) 
(4.53) 
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{Hie) - £ [B*Bl L+1 - B%^B hL+1 - S V^l+J ( 4 - 54 ) 

^ j=i+i 

Here the symbol ~ means equality up to terms vanishing in the limit L — > oo. 
In arriving at these expressions we used the fact that some terms vanish in 
the limit L — > oo: 

Lemma 4.1:: For 1 < -D < 2 the following holds in the limit L — > oo: 

(1) B 3 £ tJ^L+i - 0, (4.55) 

(2) 5 3 £ fe- 2 ^ - 0, (4.56) 

j=i+l 
L-l 

(3) 5 3 ^ &J 3 ->0. (4.57) 

j=i+l 

Proof: 
0<(1)<5 2 EU^ 2 
< (2) < fi 2 Efc+i &7 2 

Using Proposition 2.1, it is seen easily that both upper bounds go to zero for 
L — > oo. (3) is even more elementary. ■ 

There are also sums that converge to nonzero limits: 

Lemma 4.2: For 1 < D < 2 and b > 

(1) A = lim B 2 £ 67 1 S i>L+1 = i (4.58) 

ij— >oo -i-i Z 

j =1+1 

(2) h = lim B 3 £ V'^l+i = i (4-59) 

J =1+1 

For 6 = l\ and 72 vanish, because each term in the sums vanishes before 
taking the limit. 

Proof: The two statements are proven by 'summation by parts'. First 
note that for any sequences fi,gi we have 

B B 

[fi+i ~ fi\9j = - 2 - 9j]fj+i + Ib+i9b+i - f aQa- (4.60) 

j=A j=A 
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Applying this to the sum in (1) above and noting that bj 1 = Bj + i L+ i~Bj L+1 , 
we obtain 

E *] l B hL+1 = £ bf + \b? +1iL+1 - h L \ (4.61) 

j=i+l j=i+l 

From this the assertion (1) follows easily, using Lemma 4.1. 

Applying the formula ( (4.60| ) to the expression in (2) above yields, after 
some simplification, 

E bfBl L+1 = g b?B jtL¥1 - ~ E bf - \b? + \Bl hL+l . (4.62) 

j=i+l j=i+l j=i+l 

Again the assertion follows from Lemma 4.1. ■ 

Now we can add up all the contributions; the divergent terms cancel and 
we obtain 

(/) + (//) + (///) ~ {-±B 2 Bl tL+1 + (N + 1) - S(N + l)h + ANI 2 ) 

(4.63) 

For b = I\ and I2 vanish, B 2 B\ L+1 = 1 and we recover our old result 
(|4.471 ). For b > 0, however, the term B 2 B\ L+1 vanishes and the sums I\ and 
I2 converge to 1/2 and 1/3. So we obtain 

(.(.>«(i + l)) W M)= (Jy -^- 5) . (4.64) 

for b > 0. 

So we have found that in our model PT gives results depending on b.c. for 
any dimension (except for = 2)! While this is a sensible result for D > 2 
due to the occurrence of SSB, it signals a disease of PT for 1 < D < 2, where 
we have a unique Gibbs state. The right asymptotic expansion to order j3~ 2 
is obtained only with free b.c, where the thermodynamic limit is reached 
already for finite L. Free b.c. are 'DLR b.c' for our model, corresponding 
to exactly integrating out the variables outside our box. 

We now turn to the general case (s(i) • s(k))^ (a, b) for i > k. The 
computation is similar, only this time we have to split the sums into 3 parts: 
from 1 to i — 1, from i to k — 1 and from k to L or L — 1. For the sums of 
finite range we can again take the termwise limit L — > 00, which makes the 
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sums from 1 to i — 1 disappear. The sums from i to k — 1 produce, among 
others 

Using this formula is is not difficult to find the final result: 

Proposition 4.3: The one loop PT result for the spin-spin two point 
function is for free b.c. (b = 0): 



( S (z)- S (fc)) (2) (a, b) 



(N — i)(N — 3) 



2(N~ 1) 



=k-l 



-2 



(4.66) 

and represents the correct aymptotics in the thermodynamic limit. Taking 
the thermodynamic limit termwise for b > 0, we instead obtain the incorrect 
result 



( S (z)- S (fc)) (2) (a, b) 



(N — 1)(N — 5) 
8^ 



52k - 2(JV 



=fc-i 



E 



-2 



(4.67) 

A final remark concerns free b.c. at site 1 (a = 0): since the result for 
b = is strictly independent of L as well as a, and since the labeling of the 
sites can be inverted (j — > L—j), the result for a > 0, b = is identical with 
the result for a = 0, b > 0, and therefore produces the correct asymptotics. 

Comparing the results obtained with different b.c, we can sum up what 
we have found as follows: the second order PT term for the invariant two 
point function, computed by the conventional termwise thermodynamic limit 
with b.c. parameters a, b > differs from the correct second order term by 

(N — 1)(N — 2) 



5 (s(i) ■ s(k)) (2) (a,b) 



up 2 



-B 



ik 



(4.68) 



This dependence on b.c. we have found highlights a general problem of PT: 
an infrared regulator is needed, but the standard procedure of removing that 
regulator termwise yields ambiguous results. In general they are incorrect, 
and only in our simplified model we are in the lucky situation to know the 
DLR b.c. which yield the true answer. 

The case iV = 2 is special: we have not found manifest signs of a disease 
of PT. In fact for the translation invariant 0(2) model it has been proven in 



[19] that standard PT yields the correct asymptotic expansion. 
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4.5 Remarks on the magnetic field as infrared regula- 
tor 



A popular infrared regulator in the non-linear O(N) u-models is the magnetic 
field. It was used for the first detailed study of the perturbative Renormal- 
ization Group and the asymptotic freedom predicted by it . On the other 



hand it has been known for many years that the usual procedure of doing the 
perturbation expansion in the presence of a magnetic field and than remov- 
ing it termwise yields incorrect results already for a finite number of spins 
PU| and in D = 1 ||. Since our models are essentially one-dimensional, the 
same phenomenon is expected to occur also here. This is discussed in detail 
in p|. 



Here it would lead us too far afield to go into this very technical matter, 
which involves interesting methods from the theory of continued fractions. 
But it should be seen that introducing the b.c. parameters a and b essentially 
amounts to the introduction of a local magnetic field, and we have seen in 
the previous subsection that for N > 2 with a > and b > PT produces 
incorrect results, independent of the values of a and b; hence sending a, b — > 
in the end does not help. So it should be clear that one cannot hope for 
anything better with a global magnetic field. 



5 Percolation Properties 

In this section we discuss briefly some percolation properties of the model. 
Even though percolation is rather trivial in our essentially one-dimensional 
systems, we find it worthwhile to check whether the general ideas of || || |9j 
apply in this case. Not surprisingly, we find again the familiar dichotomy 
between the situation in D < 2 and D > 2. 

We are interested in the percolation properties of sets defined by the spin 
s lying in certain open connected subsets A of the sphere Sjv-i, such as the 
'polar caps' "P+ = s ■ e N > e/2 and 'equatorial strips' S e = \s ■ e N \ < e/2 
discussed in [Q, [8|, [|. For simplicity we say 'a certain subset A C SV-i 
percolates' when we mean that the set of points of our lattice {i G Z D \s(i) G 
A} percolates. 

Let us first discuss the case of no symmetry breaking (1 < D < 2): 
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Theorem 5.1: For 1 < D < 2 a subset A C S^-i whose complement is 
open and nonempty never percolates. 

Remark: If one interpretes the model as living on Z D , this can be viewed 
for D = 1 as the formation of 'islands' and for D = 2 as the 'ring formation' 
discussed in 

Proof: Assume the contrary. Consider the characteristic function XA(s{i))- 
Then by Theorem 3.2 

(Xa(s(i))) = f du(s), (5.1) 

J A 

which is a number independent of i and < 1. On the other hand, if A 
percolates, 

lim (xA(s(i))) > Um XA(s(i)) = 1 (5.2) 

by Fatou's lemma, which is a contradiction. ■ 

We now turn to the case of spontaneous symmetry breaking (D > 2). In 
this case it is to be expected that in the Gibbs state OooDk > obtained as 
the thermodynamic limit with ejv-Dirichlet b.c, there is percolation of any 
neighborhood of the 'north pole' ejy. To actually prove this requires rather 
detailed technical estimates (cf. |IIJ). Here we will limit ourselves to giving 
a simple proof of this fact for D > 3) and then show that it can be extended 
to D > 5/2. 

Theorem 5.2: For D > 5/2, in the state obtained as the thermodynamic 
limit with ejv-Dirichlet b.c, any open set A containing the point G SV-i 
percolates. 

Proof: We will show that the probabilities for the events 

a, = {s(i) i A} (5.3) 

are summable; the theorem then follows from the Borel-Cantelli Lemma, 
which states that in this case with probability 1 only finitely many of the 
events a% occur. 

We have the following 

Proposition 5.3: Let %\ be the characteristic function of the set B = 
{s(i) G Sjv-i| \s(i) — > e}. Then 

(xDocDir < « 2 " D (5-4) 
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Proof: Obviously for any n > 



tf< lfWz52i . (5.5) 



Using n = 2 and the results of Section 3 we obtain: 

i 2 \ / 5 2 \ N-l f A 2 (oo) N 



Z(Z + AT - 2) 1 , 2 _ D 
2/? 2 D D(D-2) % ' 



(5.6) 



The last expression is clearly summable over i for D > 3, so this proves 
percolation for .D > 3. 

To extend the proof to D > 5/2, we have to work a little more. Again 
by the Borel-Cantelli Lemma, the claim will follow from a sharpening of 
Proposition 5.3: 

Proposition 5.4: 

(x|) Dir < « 2(2 - D) (5.7) 



Proof: We choose n = 4 in the inequality |5.5| . Denoting s(i) • ejv = 2 we 
expand \s(i) — e^l 4 in Gegenbauer polynomials in z: 

|s(i)-e7v| 4 = 1 - 2^ 2 + 2 4 = a + a 2 C 2 2 (2) + a 4 C 4 2 (z) (5.8) 



with 

N 2 -l 



N(N + 2) ' 



(5.9) 



a 2= r ig±i) „ (5.10) 
24 

iV(iV2_4)(iV + 4)- (5 - U) 



Thus we obtain, proceeding as in Section 3 



s(i)-ejv| ) = (4>oAs - e N \ T ioo 5e N ) 

/ oo,Dir 

a + a 2 C 2 2 (l)A 2 (i, 00) + a 4 C 4 2 (l)A 2 (i,oo) (5.12) 
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Using now the asymptotics on A/(i, oo) found in Section3, we find that the 
constant terms, as well as the terms 0(i 2 ~ D ), cancel and we are left with an 
expression 0(i 2( - 2 ~ D ^) as claimed. ■ 

Remark: By developing the asymptotics of A;(z,oo) further, one could 
presumably push the percolation threshold down to 2. Alternatively, perco- 
lation for all D > 2 would follow also from the fact that the kernel T ioo (e n , s) 
decays as a Gaussian with variance (\s(i) — e 7v| 2 )ooDir- 



6 Conclusions 

The main purpose of this study of a solvable familiy of models was to test 
certain ideas of Patrascioiu and Seiler concerning the two-dimensional O(N) 
spin models. Let us review the outcome. 

The central thesis of Patrascioiu and Seiler put forward in numerous 
publications (see for instance |2|, ||, [| and references therein) is that there is 
no fundamental qualitative difference between the 'abelian' case N = 2 and 
the 'nonabelian' one N > 2. This is fully confirmed in the solvable models 
studied here: critical behavior depends only on the dimension parameter D, 
with the models becoming critical at D = 2 independent of N. 

Another point stressed by Patrascioiu and Seiler j5| is the fact that per- 
turbation theory is ambiguous in D < 2 and for N > 2; they suggest that 
the difference found in the perturbative renormalization group between the 
cases N = 2 and N > 2 is an artefact of perturbation theory. This is again 
borne out here, as discussed in Section 4. 

Finally concerning percolation properties, we find that the 'ring forma- 
tion', proposed in 0, |[ |§ for the 'soft phase' in D = 2 actually takes place 
in the models studied here. 

Due to their simplicity these models are, however, lacking certain features 
that the translation invariant O(N) models possess: due to their essentially 
one-dimensionl nature there is no difference between 'ring formation' and 
formation of 'islands' typical for the massive high temperature phase. In fact 
for D = 2 the high temperature phase has been eliminated altogether, in 
accordance with the fact the models are more ordered than the translation 
invariant ones. So this unfortunately eliminates the possibility of studying 
critical behavior in (3. 
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In our view the models studied here give support for the ideas of Pa- 
trascioiu and Seiler, but of course they cannot provide anything like a proof 
of them. 
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